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Dynamical Lie-algebraic method for the construction of 
local quantum invariants of motion in non-integrable many- 
body systems is proposed and applied to a simple but generic 
toy model, namely an infinite kicked t — V chain of spinless 
fermions. Transition from integrable via pseudo-integrable 
(intermediate) to quantum ergodic (quantum mixing) regime 
in parameter space is investigated. Dynamical phase transi- 
tion between ergodic and intermediate (neither ergodic nor 
completely integrable) regime in thermodynamic limit is pro- 
posed. Existence or non-existence of local conservation laws 
corresponds to intermediate or ergodic regime, respectively. 
The computation of time-correlation functions of typical ob- 
servables by means of local conservation laws is found fully 
consistent with direct calculations on finite systems. 

PACS numbers: 03.65. Fd, 05.30. Fk, 05.45.+b 

We investigate the existence of Local Quantum Invari- 
ants of motion (LQI) (i.e. conservation laws) in a generic 
quantum many-body system of interacting particles. For 
sufficiently strong non-linear coupling between particles 
one expects quantum mixing and quantum ergodicity 
which is incompatible with existence of LQI Q . In this 
regime (hereafter referred to as quantum ergodic), time- 
correlations between arbitrary pair of quantum observ- 
ables decay to zero. This property of quantum ergodic- 
ity is absolutely necessary for the derivation of the laws 
of normal transport within linear response theory. In 
another extreme case of completely integrable quantum 
many-body systems, an infinite number of LQI exist, 
and such systems are manifestly non-ergodic. Anoma- 
lous (ideal) transport properties of completely integrable 
quantum many body systems have been discussed in || . 
Our hypothesis is that generic non-integrable quantum 
many-body system of (locally) interacting particles may 
not be quantum ergodic, not even in thermodynamic 
limit (TL) (size — > oo and fixed density), provided that 
the system is sufficiently close to some completely inte- 
grable point in the parameter space (analogy with or- 
der to chaos transition of classical systems of few inter- 
acting particles). In a recent work 0] we conjectured 
and demonstrated the existence of intermediate, neither 
completely integrable nor ergodic, dynamical regime of 
a generic non-integrable quantum many-body system in 
TL, by directly inspecting its time evolution. It is the 
purpose of this paper to show that such an intermedi- 
ate regime may be characterized by quantum pseudo- 
integrability: existence of at least one or few LQI, which is 
a sufficient condition, using a relation proposed by Mazur 
[^| and Suzuki ||] , to prevent time correlations of certain 



observables to decay to zero. 

In Q a novel family of simple but generic many-body 
systems of locally interacting particles has been intro- 
duced smoothly interpolating between integrable and er- 
godic regime, namely kicked t- V model (KtV) of spin- 
less fermions with periodically switched nearest neigh- 
bor interaction on an infinite chain with time-dependent 
Hamiltonian 

= E [-5*( c i c i+i + h - c -) + S P {T)Vn j n j+1 ] , (1) 

3 

Cj,Cj,nj = CjCj are fermionic creation, annihila- 
tion and number operators, respectively, and 8 p (t) = 
Sm=-oo &{ T ~ m )- The time-dependent Hamiltonian (Jl|) 
can be written as H(t) = tH\ + S p (t)VH where the di- 
mensionless kinetic energy Hi and the kick potential Hq 
may be rewritten in terms of independent spin variables 
(erf , o~j) on sites j, via Wigner- Jordan transformation, 

Symmetric time evolution of KtV system for one 
period is given by explicit unitary quantum many- 
body map [h = 1), U = fexp(-if*{ 2 /2 dTH(T)) = 
exp(-itH 1 /2) aq>(-iVH Q ) exp(-itffi/2). Note that V is 
a cyclic parameter of period 2ir, and the dynamics is es- 
sentially invariant w.r.t. transformations t —* —t and 
V — » —V, so we consider only the half-strip (t, V) E 
[0, oo) x [0, tt). KtV model is completely integrable for: 
t = (Ising model), or V — (mod 2ir) (free fermion 
model), or t, V -> and t/V finite (XXZ model). 

We prefer to consider Heiseberg representation and 
write a map over the algebra il of quantum observables 
A(t) for time evolution over one period, (7 a d : A(n—^) — > 
A(n + §) = U^A(n - \)U, explicitly as 

C/ a d = exp ^^adiJi^ exp (iV&c\H ) exp ^— adiJi^ . 

where ad is the usual adjoint map on the Lie algebra il, 
a,&A:B-*[A,B]= AB - BA. Infinite-dimensional Lie 
algebra il of local quantum observables is also a Hilbert 
space equiped with the invariant bilinear form 

(A\B) = lim y^jtr^ts (2) 

which is a scalar product invariant under the adjoint 
map ad A : B -> [A,B], ((a,c\A^)B\C) = (B\(adA)C). 
The limit and the trace tr^ is defined thru finite systems 
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of increasing size L. Our aim is to check the existence 
of LQI which are the (normalizable jjj) fixed points A 
{{A\A) < oo ) of the Heissenberg map C/ a d 

U,dA = A. (3) 

However, we do not suggest to solve eq. (|3|) in the full 
algebra it which is a highly prohibitive task. Instead, 
we devise a special subalgebra of il, the Minimal In- 
variant Lie Algebra (MIL A), which is invariant to mo- 
tions generated by the kinetic or the potential part of 
the Hamiltonian and hence it is also invariant to [Tad • 
Having the two generators, Hq and H±, spanning 2-dim 
subspace 5 = {aHo + (3Hi}, we construct the basis of 
MILA & = U^Li( a ds) n s ordered by the order of locality 
as follows: We assign an observable H pb to an ordered 
pair of integers (p, b) , order p, and code b, < b < 2 P 
with p binary digits b n , b — J2n=o bn2 n , namely 

H Pi b = (ad H bp _ 1 ) (ad H bp _ 2 ) • • ■ (ad H bl )H bo . 

Since not all observables H p b upto a given maximal or- 
der q,p < q, are linearly independent we perform Gram- 
Schmit orthogonalization w.r.t. the scalar product (0) 



G = < i> c l V (Gq,c\G q . c ); G q . c ^ 0, 
I 0; G q . c = 0, 

(p,b)<(q,c) 

Gq,c — H qtC — ^ Gp tb (G Ptb \Hq tC ). (4) 
(P,b) 

The nonzero (normalized) local observables G 9jC form 
the orthonormal basis of MILA. Note that observables 
Gq. c are local operators of order q, they have been ex- 
panded in terms of spatially homogeneous finite products 
of field operators, say Z SoSl ... Sq = £\ • • • a^ q , 

where Sk G {0, +, — ,z} and cr^ = 1. The number of 
terms Z So ... Sq in expansion of G g>c was found to grow 
exponentially as ~ 2.55 9 (on average). The observables 
Z So ...s q form a convenient orthonormal Euclidean basis 
w.r.t. (||) of the Hilbert space of local spatially ho- 
mogeneous observables. Let us now consider truncated 
linear subspaces of MILA, & p = Un=o( acls ) ns ' w ^ n 
dimensions d p = dim6 p , linearly spanned by observ- 
ables Gq tC up to maximal order p, q < p. Let H Pi0 , 
denote real and symmetric (Hcrmitcan in general) ma- 
trices of linar maps ad H a on & p with images orthogo- 
nally projected back to & p . It follows from the construc- 
tion that they have (generally) a block-band structure 
where blocks correspond to observables with fixed order 
q: namely (G q . c \&d H a \G q >^ c ') ^ only if \q — q'\ = 1. 

Our program is to solve the equation ([}]) numerically in 
truncated subspaces of MILA, <S P , and check whether the 
procedure converges as p increases. This is quite feasi- 
ble since the dimensions d p = dim 6 p increase much less 
rapidly than, say, the dimensions of truncated subspaces 
of a huge Lie algebra of homogeneous observables, 



spanned by Z So ... Sp . In case of KtV the former increase 
approximately as d p « 1.68 p_1 (see table ||) while the lat- 
ter goes as ~ 4 P+1 . The truncated adjoint maps, Hp )0 , 
have nontrivial null spaces 9Tp iQ = {A G & p , [H a ,A] G 
& p +i — 6 P }, with dimensions d p a — dim^p.c which in- 
crease approximately with the same exponent oc 1.68 p . 
By means of computer algebra we managed to go as high 
as p = 14. An important observation was that a ma- 
trix 1 — exp(iitH Pi i) exp(iVH P: o) exp(i|tH P) i) possesses 
a high-dimensional null space 91 p (i, V) whose dimension 
is, for odd p, independent of parameters t, V and equal to 
the dimension of nullspace of H p .i, dimO^-i = ^22-1,1- 
Note also that for odd order of truncation p = 21 — 1, the 
elements of null space A G 9t p (t, V) are spanned by com- 
binations of odd powers of generators, i.e. (A\G2i, c ) = 0, 
which is due to time-symmetric construction of evolution 
operator U a a. However, we are seeking for LQI A G 6, 
which are normalizable, i.e. the relative norm in the 
subspace of local operators of fixed order p, defined as 
N q (A) := J2 C \{A\G q , c )\ 2 , should be a rapidly decreasing 
function of the order q, since (A| A) = N q (A) < 00. 

Only the elements of < Jlp(t, V) which are within certain 
numerical accuracy independent of the order of trunca- 
tion p are candidates for LQI. To find them we mini- 
mize a quadratic form, the relative norm at truncation 
order N P (A), i.e. we diagonalize the operator N p — 
J2 C G P ,c ( S>G PtC in the subspace $lp(t, V). In Fig.l we plot 
the relative norm N q {A m ) of the first three eigenvectors 
A m of N p with p — 13 corresponding to the smallest val- 
ues of quadratic form N p (A m ) againts the (odd) order 
q € {1,3, . . . ,p}. We give a mesh of plots for various val- 
ues of the parameters t, V. Note that N 2 i(A m ) = due 
to symmetry. In certain region of parameter space t, V 
we have found a good agreement with exponential 

N2i-i(A m ) oc exp(-s m l) (5) 

with a positive exponent s m (t, V) for the first observable 
m = 1, and in a smaller subregion of parameter space 
even for the second observable m = 2, however with a 
smaller exponent S2 ■ Positivity of s m as determined from 
exponential fit (||) for 21— 1 = 3, 5... p — 2 has been used 
as a numerical creterion of convergence (locality) of A m . 
By comparing with results for p = 11 we have checked 
that the converged observables are (almost) independent 
on the variation of the truncation order p. For the rest 
of the null space s yt p (t,V), m > 2, we found roughly a 
uniform distribution N2i-i(A m ) ~ 1 so these observables 
cannot converge as p — > 00. We conclude that this is an 
evidence for the existence of one or two LQI in certain 
region in the parameter space, roughly where t < 1.5. 
Outside this region, all exponents s m (t, V) are zero and 
none of the eigenvectors A m converges to a LQI which is 
consistent with the property of quantum ergodicity. In 
Fig. 2 we show a phase-diagram of an exponent s±(t, V) 
of LQI Ai(t, V). In table || we give explicitly the first 
few coefficients f PiC , of expansion A — Yl(p c ) fp,cG P:C of 
the one LQI for t = V = 1, and of the two LQI for 
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t = o.2,y = i §. 
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FIG. 1. The logarithmic relative norm of the first three 
candidates A m ,m = 1,2,3 for LQI, log 10 N q (A m ), is plotted 
against (odd) order q — 21 — 1, for a square mesh of parameters 
t and V. The order of truncation is p = 13. Full curves: 
m — 1, dashed curves: m = 2, and dotted curves: m = 3. 




FIG. 2. The largest inverse order-localization length si 
for the first numerical LQI Ai (p = 13) vs. parameters t, V. 
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TABLE I. Dimensions of truncated MILA 6 p , and of 
nullspaces of H p , a , for different orders of truncation p. 



Further, we use a theorem of Mazur g] generalized 
by Suzuki |J (MS) to compute canonical averages H 
of time averaged correlation functions of certain observ- 
ables. We consider dimensionless kinetic energy T = H\, 
where (T) = (1|T) = 1, with time-correlator 

1 N 

D = J im TjH^l) - < T >]IP>+ \) ( T )])- 

N—>oo 1\ z — * 
n=0 

MS equation expresses D in terms of a sum over all LQI 
D = J2\(A m \T)\ 2 /(A m \A m ). (6) 

m 

In case of complete integrability there are infinitely many 
LQI, in case of quantum ergodicity there are none and 
MS correlator D is zero, while in intermediate regime 
the sum @ (for KtV) has one or two nonzero terms. 
In Fig. 3 we show a phase diagram of the kinetic time- 
correlator D as determined from LQI (||). Note a sharp 
(phase) transition between ergodic dynamics (disordered 
phase D = 0) and intermediate dynamics (ordered phase 
D > 0) which may also be characterized by the maximal 
inverse order-localization length s\ (Fig. 2) which linearly 
decreases to zero at the transition. We have compared 
the above results on infinite KtV chains to direct calcu- 
lations on finite chains of L sites with periodic boundary 
conditions, cl+i — c\, which have a discrete quasi-energy 
spectrum r\ n and eigenstates \n), U\n) — exp(— ii] n )\n). 
For a finite system of size L, a time-correlator reads 

^ L = 2- i ^(i(n|T|n)-(T)) 2 . (7) 

n=l 

However, Dl need not necessarily converge to the proper 
time-correlator of an infinite system L = oo, since then 
the time-limit r — * oo is taken prior to the size-limit 
L — > oo Q. Nevertheless the behavior of Dl for L = 16 
shown in Fig. 4 as a function of parameters t, V is very 
similar to the infinite system MS correlator (^|) shown 
in Fig. 3. Agreement is even quantitative, except in the 
region of transition between dynamical phases (Fig. 5). 
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TABLE II. Coefficients of expansion of the two LQI 
(col.2,3) for t = 0.2, V = 1 and of the single LQI (col.4) 
for t = V = 1 in terms of basis G QtC (col.l, c is a binary code) 
up to 6th order. Truncation order is p = 13. Note that all 
the digits shown (except possibly the last one) are the same 
for truncation at p — 11. Other obsevables G q , c ,q = 1,3,5, 
are zero by construction (0). 
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Note that KtV map C/ a d is invariant under the par- 
ity operation P : cj —> C-j and MILA may exhaust |l(]] 
only the positive parity class, PA = A. Unfortunately, 
the current J — *Q3j c ] c j+i ~ h.c), which, interestingly, 
gives rise to ideal transport in intermediate [Q (inte- 
grable ||||) regime, has a negative parity PJ = — J, 
and hence zero overlap with the above LQI of MILA, 
(A m \J) = (PA m \PJ) = -(A m \J) = 0. However, analo- 
gous results have been found for the negative parity class 
of observables as well . 

The total occupation number N = J2j n j = \ J2j + 
1), is a trivial invariant of motion since (N\&) = (N\ J) = 
0. One may either consider observables over the Fock 
subspace of states with a fixed density p := (N) E [0, 1], 
or as we do here using a full trace average (g), consider 
the entire Fock space, which is in TL equivalent to half 
filling, p = §. 
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FIG. 3. MS correlator ([| D (p = 13) vs. parameters t, V. 
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FIG. 4. The finite size kinetic correlator (Q) Dl for L — 16 
vs. parameters t, V. 

Conclusions: We have found strong evidence for the ex- 
istence of non-trivial LQI, in a simple but generic quan- 
tum many-body system in TL, namely kicked t-V model 
j|. The algebraic method, which should be applicable 



to other non-integrable quantum many-body systems, is 
based on the (computerized) construction of minimal in- 
variant infinitely dimensional Lie algebra, MILA, gener- 
ated by the essential parts of Hamiltonian (in our case, 
by kinetic energy and kick potential). LQI are found 
numerically as fixed points of the adjoint map of the 
evolution operator (or of the Hamiltonian if system was 
autonomous) in MILA. Existence of LQI is found to 
be fully consistent with deviations from quantum er- 
godicity characterized by non-vanishing averaged time- 
autocorrelations D of a typical observable; here we use 
the kinetic energy. D is a suitable order parameter de- 
scribing the phase transition from the pseudo-integrable 
(intermediate) regime (D > 0) to the quantum ergodic 
regime (D = 0). 

Discussions with Prof. P. Prelovsek, and the financial 
support by the Ministry of Science and Technology of R 
Slovenia are gratefully acknowledged. 
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FIG. 5. The comparison between MS correlator (||) D of 
infinite KtV chain (thick full curve) and finite size correlator 
(|) D L for L — 16 (dotted curve) and linearly extrapolated 
to 1/L = from data for L = 16 and L = 12 (dashed curve) 
vs. parameter t, and for fixed parameter V = 1. Medium full 
curve and thin full curve denote separate contributions of the 
first and second LQI to MS correlator D m), respectively. 
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